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Abstract 

Photoproduction of beauty and charm quarks in events with at least two jets has 
been measured with the ZEUS detector at HERA using an integrated luminosity 
of 133 pb" 1 . The fractions of jets containing b and c quarks were extracted 
using the invariant mass of charged tracks associated with secondary vertices 
and the decay-length significance of these vertices. Differential cross sections 
as a function of jet transverse momentum, p£ , and pseudorapidity, if et , were 
measured. The data are compared with previous measurements and are well 
described by next-to-leading-order QCD predictions. 
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1 Introduction 



The study of beauty and charm production in ep collisions constitutes a rigorous test of 
perturbative Quantum Chromodynamics (QCD) since the heavy-quark masses provide 
a hard scale that allows perturbative calculations. At leading order, boson-gluon fusion 
(BGF), 7# — > qq with q e {b, c}, is the dominant process for heavy-quark production 
at HERA. When the negative squared four-momentum exchanged at the electron vertex, 
Q 2 , is small, the process can be treated as photoproduction, in which a quasi- real photon 
emitted by the incoming electron interacts with the proton. For heavy-quark transverse 
momenta larger than or comparable to the quark mass, next-to-leading-order (NLO) 
QCD calculations in which the massive quarks are generated in the hard sub-process [I] 
are expected to provide reliable predictions for the photoproduction cross sections. 

Beauty and charm photoproduction has been measured using several different methods 
by both the ZEUS and HI collaborations. In most of the previous measurements of 
beauty photoproduction at HERA, the cross section was determined using semileptonic 
decays into muons [2H5] or electrons [HIE]- In the muon analyses, the fraction of leptons 
originating from beauty was determined by using the large transverse momentum of the 
muon relative to the axis of the associated jet, p^ 1 , and/or exploiting the impact parameter 
of the muons. In the more recent electron analysis j6], several variables, sensitive to both 
electron identification as well as to semileptonic decays, were combined in a likelihood- 
ratio test function in order to extract the beauty and charm content. The HI collaboration 
has published an inclusive measurement of beauty- and charm-quark photoproduction 
using a method based on the impact parameter of tracks to the primary vertex [5]. The 
other published charm or beauty photoproduction measurements [9HT5] used either meson 
tags or a combination of lepton and meson tags. In all of the above analyses reasonable 
agreement between the measurement and the theory prediction was found. 

The aim of this measurement is to test perturbative QCD with high precision. For this 
purpose, the long lifetimes of the weakly decaying b and c hadrons as well as their large 
masses were exploited. The measurement relies on the reconstruction of decay vertices 
with the ZEUS silicon microvertex detector (MVD) [16J. Two discriminating variables 
were used: the significance of the reconstructed decay length and the invariant mass of the 
charged tracks associated with the decay vertex (secondary vertex). The measurement 
was kept fully inclusive, leading to a reduced uncertainty due to branching fractions and 
a substantial increase in statistics compared to exclusive analyses. The high statistics 
also allowed the kinematic region of the measurement to be extended to high values of 
the transverse jet momentum, . 
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2 Experimental set-up 



The analysis was performed with data corresponding to an integrated luminosity of 
133 pb _1 which were taken during 2005. Electrons at an energy of E e = 27.5 GeV collided 
with protons at E p = 920 GeV, yielding a centre-of-mass energy of 318 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [T7j. A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central track- 
ing detector (CTD) [18] and the microvertex detector (MVD) [16]. These components 
operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The 
CTD consisted of 72 cylindrical drift-chamber layers, organised in nine superlayers cover- 
ing the polar-angle^ region 15° < 9 < 164°. The MVD silicon tracker consisted of a barrel 
(BMVD) and a forward (FMVD) section. The BMVD contained three layers and pro- 
vided polar-angle coverage for tracks from 30° to 150°. The four-layer FMVD extended the 
polar- angle coverage in the forward region to 7°. After alignment, the single-hit resolution 
of the MVD was 24 /im. The transverse distance of closest approach (DCA) to the nominal 
vertex in X—Y was measured to have a resolution, averaged over the azimuthal angle, of 
(46©122/pr) /im, with p? in GeV. For CTD-MVD tracks that pass through all nine CTD 
superlayers, the momentum resolution was ct(pt)/pt = 0.0029pT © 0.0081 © 0.0012/pt, 
with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [Tj5] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a(E)/E = 0.18/ \/E for electrons and 
a(E)/E = 0.35/ VE for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — > cyp by a luminosity 
detector which consisted of independent lead-scintillator calorimeter [20] arid magnetic 
spectrometer [21] systems. The fractional systematic uncertainty on the measured lumi- 
nosity was 1.8 %. 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the centre of the CTD. The pseudorapidity is defined 
as r\ = — In (tan |), where the polar angle, 9, is measured with respect to the proton beam direction. 
The azimuthal angle, <f>, is measured with respect to the X axis. 
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3 Monte Carlo simulation 



Monte Carlo (MC) samples of beauty, charm and light-flavour events generated with 
Pythia 6.2 [22] were used to evaluate the detector acceptance and to provide the predic- 
tions of the signal and background distributions. 

The production of bb and cc pairs was simulated following the standard Pythia prescrip- 
tion, using leading-order matrix elements combined with parton showering. The following 
subprocesses [23] were generated: 

• direct and resolved photoproduction with leading-order massive matrix elements. In 
the direct-photon process, the quasi-real photon enters directly in the hard interaction, 
while in the resolved-photon process, the photon acts as a source of light partons which 
take part in the hard interaction. The fe-quark and c-quark masses were set to 4.75 GeV 
and 1.5 GeV, respectively; 

• 6-quark and c-quark excitation, i.e. the contribution to the leading-order massless 
matrix elements of b and c quarks from initial-state photon or gluon splitting. 

The light-quark predictions were taken from a simulation of both direct and non-direct 
inclusive photoproduction with leading-order matrix elements in the massless scheme. 
This sample also includes final-state gluon splitting into bb and cc pairs, which is treated 
as part of the signal. 

The CTEQ4L [24] and CTEQ5L [25] proton parton distribution functions (PDFs) were 
used for the light-flavour and heavy-flavour samples, respectively. The GRV-G LO [26] 
photon PDF was used for all samples. 

The lifetimes of the B ± , B° and B s mesons were corrected from the default Pythia values 
to reflect the world averages [27]. 

The generated events were passed through a full simulation of the ZEUS detector based 
on Geant 3.21 [28]. The final MC events had to fulfil the same trigger requirements and 
pass the same reconstruction programme as the data. 

4 Data selection and event reconstruction 

A three-level trigger system was used to select events online [T71[29l[30] . At the third level, 
jets were reconstructed using the energies and positions in the CAL. Events with at least 
two jets with transverse momentum in excess of 4.5 GeV within \r}\ < 2.5 were selected. 

The tracking efficiency at the first-level trigger (FLT) as well as the efficiency of the 
dijet third-level trigger (TLT) were lowered in the detector simulation such that they 
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reproduced the efficiencies as measured in the data. The trigger efficiencies were ~ 86 % 
for the FLT and 76-100 % for the TLT, depending on the transverse momentum of the 
jets, with an average of about 90 %. The average corrections amounted to ~ 7.7 % for the 
FLT and « 3.7% for the TLT. 

The hadronic system was reconstructed from energy-flow objects (EFOs) [31] combining 
track and calorimeter information, corrected for energy loss in the dead material. Each 
EFO, i, was assigned a reconstructed four-momentum (px,PyiPzi assuming the pion 
mass. Jets were reconstructed from EFOs using a kx clustering algorithm [32] in the 
longitudinally invariant mode [33]. The ^-recombination scheme, which produces massive 
jets whose four-momenta are the sum of the four-momenta of the clustered objects, was 
used. 

At least two jets with \rf ct \ < 2.5 and pij, > 7(6) GeV for the highest (second highest) 
energetic jet were required. Only events with a well reconstructed primary vertex with 
|-Zvtx| < 30 cm were selected. 

In order to remove background from deep inelastic scattering (DIS), events were rejected in 
which a scattered-electron candidate was found in the calorimeter with energy E' e > 5 GeV 
and y e < 0.9, with y e — 1 — (1 — cos^), where 6' e is the polar angle of the outgoing 
electron. The event inelasticity, y, was reconstructed from the hadronic final state using 
the Jacquet-Blondel method [31] with ?/jb = — p l z )/2E e , where the sum runs over 

all the EFOs. A cut 0.2 < yjB < 0.8 was used to remove residual DIS events and non-ep 
interactions. These requirements correspond to an effective cut of Q 2 < 1 GeV 2 with a 
median of Q 2 fa 10~ 3 GeV 2 , as estimated from simulations. 

In order to reconstruct secondary vertices related to b- and c-hadron decays, tracks were 
selected if: 

• p T > 0.5 GeV; 

• the number of superlayers in the CTD > 3; 

• the total number of hits^l in the MVD > 4. 

The tracks were associated with one of the two highest energetic jets if they fulfilled AR = 

or more of such tracks were associated with the 
selected jet, a candidate vertex was fitted from the selected tracks using a deterministic 
annealing filter [35]. This fit provided the vertex position including its error matrix as 
well as the invariant mass, m vtx , of the charged tracks associated with the reconstructed 
vertex. Vertices with x 2 /ndf < 6, a distance from the interaction point within 1 cm in the 
X—Y plane and ±30 cm in the Z direction, and 0.8 < m vtx < 7.5 GeV were retained for 
further analysis. 

2 Each MVD layer provided two coordinate measurements. 
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Only those secondary vertices that were associated with one of the two jets with the 
highest were considered, since these jets were most likely to correspond to heavy- 
quark jets. The associated jet was required to be reconstructed within the central part of 
the detector with —1.6 < rf et < 1.4. 

5 Extraction of the heavy-flavour cross sections 

Using the secondary vertices associated with jets, the decay length, d, was defined as the 
distance in X—Y between the secondary vertex and the interaction point|f|, projected onto 
the jet axis in the X—Y plane. 

The decay-length significance, S, was defined as d/5d, where 5d is the uncertainty on d. 
The sign of the decay length was assigned using the axis of the jet to which the vertex is 
associated: if the decay-length vector was in the same hemisphere as the jet axis, a positive 
sign was assigned to it; otherwise the sign of the decay length was negative. Negative 
decay lengths, which originate from secondary vertices reconstructed on the wrong side 
of the interaction point with respect to the direction of the associated jets, are unphysical 
and caused by detector resolution effects. A small correction [23] to the MC decay-length 
distribution was applied in order to reproduce the negative decay-length data: 5 % of 
the tracks in the central region were smeared and an additional smearing was applied to 
tracks in the tails of the decay-length distribution. 

The shape of the decay-length significance distribution together with the secondary- vertex 
mass distribution, m vtx , is used to extract the beauty and charm content. The invariant 
mass of the tracks fitted to the secondary vertex provides a distinguishing variable for jets 
from b and c quarks, reflecting the different masses of the b and c hadrons. Figure [T] shows 
the decay-length significance, S, divided into the three mass bins 0.8 < m vtx < 1.4 GeV, 
1.4 < m vtx < 2 GeV and 2 < m vtx < 7.5 GeV. The MC simulation provides a good 
description of the data in all three bins and an almost pure beauty region can be obtained 
at high significances in the bin 2 < m vtx < 7.5 GeV. 

In order to minimise the effect of the light-flavour contribution, the contents of the neg- 
ative bins of the significance distribution, N(S~), were subtracted from the contents of 
the corresponding positive bins, N(S + ), yielding a subtracted decay-length significance 
distribution. An additional advantage of this subtraction is that symmetric systematic 
effects, which might arise from discrepancies between the data and the MC, are removed. 

3 In the X—Y plane, the interaction point was defined as the centre of the beam ellipse, determined 
using the average primary vertex position for groups of a few thousand events, taking into account the 
difference in angle between the beam direction and the Z direction. The Z coordinate was taken as 
the Z position of the primary vertex of the event. 
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In order further to reduce the uncertainty due to remaining differences between data and 
MC in the core region of the significance distribution, a cut of | | > 3 was applied. As a 
consistency check this cut was varied in order to estimate the uncertainty due to the MC 
modelling of the low \S\ region; effects smaller than 1 % on the beauty results and 3 % on 
the charm results were found. 

After all selection cuts, a sample of 70 433 jets with associated secondary vertices re- 
mained. 

Figure [2] shows the data and MC distributions of p!£ , rf ct , ?rivtx, the secondary vertex 
track multiplicity, n tr k, and x 2 /ndi of the secondary vertices. All distributions are shown 
after all selection cuts, except for the xV n df distribution, where the x 2 /ndi cut has not 
been applied yet. Also shown in Fig. [2] is the fraction of the total hadronic E — pz carried 
by the two highest-p T jets, 

E-pz • 

weighted by the number of jets with associated secondary vertices in the event. This 
distribution is sensitive to the fraction of direct and non-direct photoproduction contri- 
butions. The MC provides an adequate description of the data for all variables except 
rf et ; the effect of this discrepancy on the results is discussed in Section 

The beauty and charm contributions were extracted using a least-squares fit [23,36j to the 
subtracted distributions in the three mass bins. The MC beauty, charm and light-flavour 
contributions, normalised to the data luminosity, were scaled by the factors k c and 
k\{, respectively, to give the best fit to the observed subtracted distributions. The overall 
MC normalisation was constrained by requiring it to be consistent with the normalisation 
of the data in the significance distribution with \S\ > 3 and 0.8 < m vtx < 7.5 GeV. 
The subtracted and fitted distributions for the three mass bins are shown in Fig. [3J 
The contribution of the light flavours was substantially reduced through the subtraction. 
After the subtraction, good agreement was also observed between the data and the MC 
simulation. The fit procedure was repeated in different bins of p!£ and if et to obtain the 
differential cross-sections da/dp^ and da/dr]^ ct . 

In order to check the quality of the data description by the MC, subtracted distributions 
of , if et , mvtx, the secondary-vertex track multiplicity, nt r k, and \S\ are shown in Fig. |4] 
after beauty enrichment (2 < m vtx < 7.5 GeV and \S\ > 8) and in Fig. [5] after charm 
enrichment (0.8 < m vtx < 2 GeV). 

The total visible cross section for inclusive heavy-quark jet production, a q , with q e {b, c} 
is given by 

a q = -jl— . 

•A-q ' *-"Data 
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Here, £oata denotes the integrated luminosity, A q is the acceptance and jVg ec ' Data the 
number of reconstructed heavy-quark jets in data, which was determined from the fit 
using 

jyrec,Data j\j-rec,WlC 

with jVg ec ' MC being the number of reconstructed events in a MC sample with the same 
integrated luminosity as the data. k q denotes the heavy-quark scaling factor obtained 
from the fit. Defining the acceptance as 

rec,MC 



the cross section can be written as 



A = —1 

^1 jytrue.HL ' 



kg ■ TV*™ 6 ' 111 " 



£Data 

Here, A r * rue ' HL denotes the number of generated heavy-quark jets at hadron level (HL). 
Hadron-level jets were obtained by running the kr clustering algorithm in the same mode 
as for the data with the E'-recombination scheme. The algorithm was run on all final-state 
MC particles before the decay of the weakly decaying b or c hadrons. True b or c jets 
were then defined as all hadron-level jets containing a b or c hadron. Signatures with b 
or c hadrons resulting from final-state gluon splitting (g — > qq) were also included in the 
respective signal, independent of the quark flavours in the hard subprocess. The contri- 
bution of gluon splitting to the beauty signal amounted to m 2 %, while the contribution 
to the charm signal was ~ 10 %. 

The single-differential heavy-quark jet production cross section as a function of a given 
variable, v, is defined accordingly: 

da* _ k q ■ ATf ue < HL 
dv ~ £ Data • Av ' 

where Av is the width of the bin. 



6 Systematic uncertainties 

Systematic uncertainties were evaluated by appropriate variations of the MC simulation. 
The fit of the subtracted decay-length significance in m vtx bins was repeated and the cross 
sections were recalculated. The uncertainties on the total cross sections determined for 
each source are summarised in Table [TJ The following sources of experimental systematic 
uncertainties were identified [23] : 
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1. the systematic uncertainties associated with the TLT and FLT trigger efficiency correc- 
tions (see Section H]) were determined by varying each correction within its estimated 
uncertainty; 

2. the calorimetric part of the jet energy was varied by ±3 %; 

3. the track-finding inefficiency in the data with respect to the MC was estimated to be 
at most 2 %. The overall uncertainty due to this tracking inefficiency was determined 
by randomly rejecting 2 % of all tracks in the MC and repeating the secondary vertex 
finding and all subsequent analysis steps; 

4. the uncertainty due to the smearing procedure was estimated by varying the fraction 
of secondary vertices for which the decay length was smeared by ±2 %. For variations 
of the fraction in this range the agreement between data and MC remained reasonable; 

5. the uncertainty due to the asymmetry of the light-flavour content of the sample was 
evaluated by varying fcif by ±11 %. The size of the variation was estimated from the 
uncertainty on the light-flavour fraction as determined by a fit to the subtracted decay- 
length significance distribution, where the overall normalisation constraint using the 
unsubtracted distribution was not applied; 

6. the MC distributions for both light and heavy flavours were reweighted in rf et and 

to account for the differences between data and MC (see Fig. [2]). A reweighting 
of only the light-flavour content was also investigated. No significant change of the 
cross sections was observed and therefore no additional systematic uncertainty was 
assigned; 

7. the various D mesons have different lifetimes and decay modes. In order to account 
for the uncertainty of the different fragmentation fractions, the D + /D° and D + /Df 
ratios were varied by ±10% while keeping the total number of c hadrons constant; 

8. the charm fragmentation function was varied by weighting all events according to 

Z = r 

\-C J ± Je-quark jet 

calculated in the string rest-frame [22] such that the change in the mean value of z 
corresponded to the measured uncertainty [37]; 

9. the beauty fragmentation function was varied in analogy to the charm case using a 
variation of z corresponding to a variation of the Peterson fragmentation parameter, 
e b , of ±0.0015 [381132!; 

10. a 1.8 % overall normalisation uncertainty was associated with the luminosity measure- 
ment. It was included in the systematic error on the total cross sections, but not in 
those of the differential cross sections. 
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The same variations were applied to each bin for the differential cross sections. The total 
systematic uncertainty was obtained by adding the above contributions in quadrature. In 
the case of beauty, the dominant effects arise from the variation of the trigger-efficiency 
corrections, the track- finding efficiency and the reweighting as a function of p J ^\ For 
charm, the variation of the trigger-efficiency corrections as well as the energy-scale varia- 
tion contribute most to the total systematic uncertainty. 

As an additional consistency check, the contributions of direct and non-direct photon 
processes were investigated by reweighting the x J 7 et distributions based on MC and data 
comparisons of the b- and c-enriched samples. The effect on the cross sections was smaller 
than that due to the reweighting of the p!^ and rf et distributions and so a further contribu- 
tion was not added to the systematic uncertainties. A reweighting of the m vtx distribution 
was also done in order to account for residual differences between the data and the MC. 
Its effect on the cross sections was found to be negligible. 



7 Theoretical predictions and uncertainties 

The measured total and differential cross sections were compared to NLO QCD predictions 
calculated with the FMNR programme [40]. This calculation is based on the the fixed- 
flavour-number scheme, using three light flavours for the charm predictions and four for 
beauty. The PDFs were taken from CTEQ6.6 [41] for the proton and GRV-G HO [26] for 
the photon. The heavy-quark masses (pole masses) were set to m& = 4.75 GeV and m c = 
1.5 GeV. The QCD scale, A^ D , was set to 0.226 GeV. The renormalisation scale, /zr, and 

the factorisation scale, fip, were chosen to be equal and set to = /ip = \ Jvt + m b(c) > 
where Pt is the average transverse momentum of the heavy quarks. In order to ease the 
comparison with previous analyses, the theoretical predictions were also made using the 
CTEQ5M [25] proton PDFs. 

For the systematic uncertainty on the theoretical prediction, the masses and scales were 
varied separately and the effects of both variations were added in quadrature. The masses 
were varied using the values = 4.5 and 5.0 GeV, m c = 1.3 and 1.7 GeV; the scales were 
varied using fiR = /ip = \ \Jf>T + m t(c) anc ^ Jpt + m 6(c) ■ ^he resu lting uncertainties on 
the NLO QCD predictions for the total cross sections are +22% and —15% for beauty 
and +42 % and —21 % for charm. 

Parton-level jets were found by applying the kx clustering algorithm to the generated 
partonic final state in the same mode as for the hadron level in the MC (see Section [5]). 
The NLO QCD predictions for parton-level jets were corrected for hadronisation effects. A 
bin-by-bin procedure was used whereby da = do"NLO ■ Chad, and d<TNLO is the cross section 
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for partons in the final state of the NLO calculation. The hadronisat ion- correct ion factors, 
Chad) were obtained from the ratio of the hadron-level to the parton-level MC jet cross 
section, where the parton level is defined as being the result of the parton-showering stage 
of the simulation. The correction factors are given in Tables [2] and El their uncertainty 
was negligible in comparison to the other theoretical uncertainties [2]. 



8 Results 

The total and single-differential beauty- and charm-jet cross sections were measured for 
the processes 

e~p ->■ e~ b(b) X 
e~p —> e~ c(c) X 

in events with 

Q 2 < 1 GeV 2 , 0.2 < y < 0.8, p^ m > 7(6) GeV, -2.5 < ?? jctl(2) < 2.5 . 

Here rf et1 ^ and p^ 11 ^ refer, respectively, to the pseudorapidities and the transverse mo- 
menta of the two jets in the event with the largest transverse momentum within the range 
\f]i ct \ < 2.5. The cross sections are measured for those jets among these two satisfying 

-1.6 < 7/"- jct < 1.4, 

with q G {b, c}. 

The total beauty- and charm-jet production cross sections were measured as 

a^ is = 682 ± 21(stat.)+ |(syst.) pb, 
a™ = 5780 ± 120(stat.)i^(syst.) pb. 

The errors given correspond to the statistical uncertainties and the total systematic un- 
certainties including the errors due to the uncertainty in the luminosity measurement. 
The measurements were compared to NLO QCD predictions calculated with the FMNR 
programme using the specifications given in Section [7) 

^ NL °®^ ad = 7401™ Pb, 
^ NLO ®^ ad = 6000^pb. 

Hadronisation corrections of C*£ ad = 0.84 and C^ ad = 0.83 were applied to the NLO QCD 
predictions. Good agreement between the measured cross sections and the NLO QCD 
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predictions is observed. Replacing CTEQ6.6 by CTEQ5M as proton PDF reduces the 
theory predictions by ~ 5 %. 

The beauty and charm cross sections as a function of p 3 ^, and 77 jet are given in Tables [2] 
and El respectively, and are shown in Fig. [6j The measurements are compared to the 
NLO QCD predictions and to the Pythia MC scaled (see Section [5]) by a factor of 
1.11 for beauty and 1.35 for charm, as obtained from the inclusive fit. The NLO QCD 
predictions are in good agreement with the data and the scaled Pythia MC describes 
the distributions well. 

In Fig. [7] the fe-jet cross section, da/drf et , is compared to a previously published analysis 
(12] using semileptonic decays into muons in dijet events. Both measurements agree well. 
The improved precision of this analysis can be clearly seen. While a direct comparison 
with a previous HI measurement using a similar approach [8] is not possible, as the cross- 
section definitions are different, the relative errors on the measurements in this paper are 
approximately a factor 3 (2) smaller for beauty (charm). 

In order to enable direct comparisons with other ZEUS measurements given at the b- 
quark level [SHIlElIZlEj, the NLO QCD prediction corrected for hadronisation was used 
to extrapolate the dijet cross sections to inclusive 6-quark cross sections: 

/ \ vis 

da UfQ /da\ NLO 
dp b T /^\ NL ° \dp b T ) ' 

For the previous measurements, the extrapolations have been updated using the CTEQ6.6 
proton PDFs. In Fig. [Sj, the 6-quark differential cross sections as a function of the quark 
transverse momentum, da(ep — > bX)/dp b r , are shown for 6-quark pseudorapidity in the 
laboratory frame, |^| < 2, for Q 2 < 1 GeV 2 and 0.2 < y < 0.8. The b quark was not taken 
into account in the definition of the 6-quark cross section. The measurement presented 
here extends the kinematic region to higher p h T values than previous measurements and 
represents the most precise measurement of 6-quark photoproduction at HERA. Good 
agreement with the NLO QCD prediction is observed for many independent ZEUS mea- 
surements, giving a consistent picture of 6-quark photoproduction over a wide range of 

Pi 

The corresponding c-quark cross sections were also calculated and are shown in Fig. 
Due to the lower mass of the charm quark, its momentum is more affected by gluon 
radiation. Hence the corresponding cross section is shown as a function of the parton- 
level jet momentum (calculated as in Section [7j) rather than that of the quark. Here the 
cross sections have been extrapolated to the region \rj c .j C t\ < 1-5, as it corresponded better 
to the measurements. 
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The c-quark jet cross sections are consistent with previous ZEUS measurements [TtfTU] 
and are in good agreement with the NLO QCD prediction. 

9 Conclusions 

Inclusive beauty- and charm-jet cross sections in photoproduction at HERA have been 
presented, exploiting the long lifetimes and large masses of b and c hadrons. Compared to 
previous measurements of specific decay chains, this analysis has substantially increased 
statistics and a reduced dependence on the branching fractions. The background from 
light-quark jets was suppressed by using the subtracted decay-length significance distri- 
bution of secondary vertices. 

The visible cross sections as well as differential cross sections as a function of p 3 ^ 1 and rf et 
have been compared with NLO QCD calculations. Good agreement is observed. 

In order to be able to compare these cross sections with others, they have been extrapo- 
lated to the region j^l < 2 (|^ c -jet| < 1-5) using the NLO QCD predictions. Cross sections 
as a function of the transverse momentum of the b quark and of the c-quark jet have been 
determined and compared with previous ZEUS measurements. The measurements agree 
with each other and give a consistent picture of heavy-quark photoproduction over a wide 
kinematic range. 

The charm cross sections presented in this paper are more precise than previous measure- 
ments made by the ZEUS collaboration and have similar accuracy as measurements made 
by HI. The beauty cross sections represent the most precise measurements of 6-quark 
photoproduction made at HERA. 

Acknowledgements 

We appreciate the contributions to the construction and maintenance of the ZEUS de- 
tector of many people who are not listed as authors. The HERA machine group and the 
DESY computing staff are especially acknowledged for their success in providing excel- 
lent operation of the collider and the data-analysis environment. We thank the DESY 
directorate for their strong support and encouragement. 



12 



References 

[1] S. Frixione, P. Nason and G. Ridolfi, Nucl. Phys. B 454, 3 (1995); 
S. Frixione et al, Phys. Lett. B 348, 633 (1995). 

[2] ZEUS Collab., S. Chekanov et al, Phys. Rev. D 70, 12008 (2004). 

[3] ZEUS Collab., S. Chekanov et al., JHEP 02, 032 (2009). 

[4] ZEUS Collab., S. Chekanov et al, JHEP 04, 133 (2009). 

[5] HI Collab., A. Aktas et al, Eur. Phys. J. C 41, 453 (2005). 

[6] ZEUS Collab., J. Breitweg et al, Eur. Phys. J. C 18, 625 (2001). 

[7] ZEUS Collab., S. Chekanov et al, Phys. Rev. D 78, 072001 (2008). 

[8] HI Collab., A. Aktas et al, Eur. Phys. J. C 47, 597 (2006). 

[9] ZEUS Collab., J. Breitweg et al, Eur. Phys. J. C 6, 67 (1999). 
[10] ZEUS Collab., S. Chekanov et al, Nucl. Phys. B 729, 492 (2005). 
[11] ZEUS Collab., S. Chekanov et al, Eur. Phys. J. C 44, 351 (2005). 
[12] ZEUS Collab., S. Chekanov et al, Eur. Phys. J. C 50, 299 (2007). 
[13] ZEUS Collab., S. Chekanov et al, Phys. Lett. B 565, 87 (2003). 
[14] HI Collab., A. Aktas et al., Eur. Phys. J. C 50, 251 (2006). 
[15] HI Collab., A. Aktas et al, Phys. Lett. B 621, 56 (2005). 
[16] A. Polini et al., Nucl. Instr. and Meth. A 581, 656 (2007). 

[17] ZEUS Collab., U. Holm (ed.), The ZEUS Detector. Status Report (unpublished), 
DESY (1993), available on http://www-zeus.desy.de/bluebook/bluebook.html. 

[18] N. Harnew et al., Nucl. Instr. and Meth. A 279, 290 (1989); 
B. Foster et al., Nucl. Phys. Proc. Suppl. B 32, 181 (1993); 
B. Foster et al, Nucl. Instr. and Meth. A 338, 254 (1994). 

[19] M. Derrick et al., Nucl. Instr. and Meth. A 309, 77 (1991); 
A. Andresen et al., Nucl. Instr. and Meth. A 309, 101 (1991); 
A. Caldwell et al., Nucl. Instr. and Meth. A 321, 356 (1992); 
A. Bernstein et al, Nucl. Instr. and Meth. A 336, 23 (1993). 

[20] J. Andruszkow et al, Preprint DESY-92-066, DESY, 1992; 
ZEUS Collab., M. Derrick et al, Z. Phys. C 63, 391 (1994); 
J. Andruszkow et al., Acta Phys. Pol. B 32, 2025 (2001). 

[21] M. Helbich et al., Nucl. Instr. and Meth. A 565, 572 (2006). 



13 



[22] T. Sjostrand et al., Comp. Phys. Comm. 135, 238 (2001); 

E. Norrbin and T. Sjostrand, Eur. Phys. J. C 17, 137 (2000); 

T. Sjostrand, L. Lonnblad, and S. Mrenna, Preprint hep-ph/0 108264, 2001. 

[23] V. Schonberg, Ph.D. Thesis, Universitat Bonn, Bonn, Germany, Report 
BONN-IR-2010-05, 2010, available on 
http : //hss .ulb.uni-bonn.de/diss_online. 

[24] H.L. Lai et al., Phys. Rev. D 55, 1280 (1997). 

[25] H.L. Lai et al., Eur. Phys. J. C 12, 375 (2000). 

[26] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 46, 1973 (1992); 
M. Gliick, E. Reya and A. Vogt, Phys. Rev. D 45, 3986 (1992). 

[27] Particle Data Group, C. Amsler et al, Phys. Lett. B 667, 1 (2008). 

[28] R. Brun et al, geant3, Technical Report CERN-DD/EE/84-1, CERN, 1987. 

[29] PD. Allfrey et al, Nucl. Instr. and Meth. A 580, 1257 (2007). 

[30] W.H. Smith, K. Tokushuku and L.W. Wiggers, Proc. Computing in High-Energy 
Physics (CHEP), Annecy, France, C. Verkerk and W. Wojcik (eds.), p. 222. 
CERN, Geneva, Switzerland (1992). Also in preprint DESY 92-150B. 

[31] G.M. Briskin, Ph.D. Thesis, Tel Aviv University, Report DESY-THESIS 1998-036, 
1998. 

[32] S. Catani et al., Nucl. Phys. B 406, 187 (1993). 

[33] S.D. Ellis and D.E. Soper, Phys. Rev. D 48, 3160 (1993). 

[34] F. Jacquet and A. Blondel, Proceedings of the Study for an ep Facility for Europe, 
U. Amaldi (ed.), p. 391. Hamburg, Germany (1979). Also in preprint DESY 79/48. 

[35] K. Rose, E. Gurewitz and G.C. Fox, Phys. Rev. Lett. 65, 945 (1990); 
K. Rose, Proceedings of the IEEE, Vol. 86, pp. 2210-2239. (1998); 

F. Didierjean, G. Duchene and A. Lopez-Martens, Nucl. Instr. and Meth. 615, 188 
(2010). 

[36] A. Yagiies, Ph.D. Thesis, Humboldt Universitat zu Berlin, Berlin, Germany, 2008, 
available on http : / / edoc . hu-berlin . de. 

[37] ZEUS Collab., S. Chekanov et al., JHEP 04, 082 (2009). 

[38] C. Peterson et al, Phys. Rev. D 27, 105 (1983). 

[39] P. Nason and C. Oleari, Nucl. Phys. B 565, 245 (2000). 

[40] S. Frixione et al., Nucl. Phys. B 412, 225 (1994). 

[41] J. Pumplin et al., JHEP 07, 012 (2002). 



14 



[42] ZEUS Collab., S. Chekanov et al., Phys. Rev. D 70, 12008 (2004). 



15 



Source 


Beauty / Charm 
(%) 


la) TLT trigger efficiency 


±0.8 / ±2.0 


lb) FLT trigger efficiency 


+4.1 / +4.0 

q a / Q V 


2) CAL hadronic energy scale 


±0.6 / ±4.3 


3) Track-finding uncertainty 


±5.9 / ±1.0 


4) Decay-length smearing 


±1.0 / ±0.7 


5) Light-flavour asymmetry 


±0.2 / ±0.7 


6a) rf et reweighting 


-1.2 / -1.0 


6b) p!^ reweighting 


-5.5 / -1.1 


7a) D ± /D° ratio 


+0 / +0.6 
-1.3 / -1.8 


7b) D ± /Df ratio 


+0 / +0.1 
-1.2 / -1.3 


8) Charm fragmentation 


+0.3 / +1.2 
-0.3 / -1.3 


9) Beauty fragmentation 


+1.8 / +0.1 
-2.1 / -0.1 


10) Luminosity measurement 


±1.8 / ±1.8 


Total 


+7.8 / +6.7 
-7.7 / -7.0 



Table 1: Systematic uncertainties on the total beauty- and charm-jet cross sections. 
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jet 

Ft 


dab/dp^ 


do^o/dp? ® CL d 


rib 
°had 


(GeV) 


(pb/GeV) 


(pb/GeV) 




6 


11 


95.6 ±4.9 ito 


109 


0.83 


11 


16 


24.8 ±1.2 t\i 


29.1 1 I;? 


0.89 


16 


21 


6.02 ± 0.49lg;p 


7 1 +2.0 
' - 1 - 1.2 


0.92 


21 


27 


0.93 ± 0.221°;^ 


1 07+ 0.54 
0.34 


0.95 


27 


35 


0.30 ± 0.121°;^ 


46+ °' 13 

u -^°- 0.08 


1.05 







dcr fe /d?7 jet 
(Pb) 


d(J NLO /dr?j et ,3 
(pb) 


°had 


-1.6 : - 


-1.1 


57 ± 22t l \ 


79+22 

' z -13 


0.70 


-1.1 : - 


-0.8 


121±2llJ{j 


18 9+50 


0.78 


-0.8 : - 


-0.5 


214 ± 22±?| 




0.79 


-0.5 : - 


-0.2 


233 ± 21±|f 


307ljjjj 


0.79 


-0.2 : 


0.1 


264 ± 22+| 


342+^ 


0.81 


0.1 : 


0.5 


316 ±21+^ 


34611? 


0.86 


0.5 : 


1.4 


288 ± 15±|g 


2651^ 


0.93 



Table 2: Summary table of differential beauty-jet photoproduction cross sections, 
as defined in Section^ The measurements are given together with their statistical 
and systematic uncertainties. The NLO QCD predictions using CTEQ6.6 and their 
uncertainty are also listed. The last column gives the hadronisation correction 
factors, C b had . 
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iet 

Ft 


daj d^T 


dof*/d# (8) Q ad 


nc 
°had 


(GeV) 


(pb/ GeV) 


(pb/ GeV) 




6 


11 


906 ±24 tlo 


967 + 380 


0.82 


11 


16 


194 ± 7 tjo 


192 1 Jj 


0.90 


16 


21 


39.1 ± 3.3+ H 


38.5 1 1 5 


0.92 


21 


27 


10.5 ± 2.1+ to 


8.9 + ^ 


0.90 


27 


35 


0.9 ± 0.71 oi 


1-96+ 8^3 


0.91 









dcr c / d?7 je1: 
(Pb) 


d(T NLO /d7?je t g cu 
(pb) 


°had 


-1 


6 : - 


-1.1 


499 ± 791 ^ 


8251 g° 


0.71 


-1 


1 : - 


-0.8 


1380 ± 1101}}° 


1 OQQ+ 700 
iyo °- 400 


0.79 


-0 


8 : - 


-0.5 


2090 ± 1201}^ 


25661 £g 


0.80 


-0 


5 : - 


-0.2 


2460 ± 1301}^ 


29481^™ 


0.80 


-0 


2 : 


0.1 


2920 ± 13011$ 


29751^™ 


0.83 





1 : 


0.5 


2600 ± HOl^g 


26021^™ 


0.87 





5 : 


1.4 


2040 ± 91±^o 


15791 $ 


0.89 



Table 3: Summary table of differential charm-jet photoproduction cross sections, 
as defined in Section The measurements are given together with their statistical 
and systematic uncertainties. The NLO QCD predictions using CTEQ6.6 and their 
uncertainty are also listed. The last column gives the hadronisation correction 
factors, C c had . 
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ZEUS 



0.8 < m„ tx < 1 .4 GeV 1 .4 < m„„ < 2 GeV 




S 



Figure 1: Distributions of decay-length significance, S, for (a) 0.8 < m vtx < 
1.4 GeV, (b) 1.4 < m vtx < 2 GeV and (c) 2 < m vtx < 7.5 GeV. The data are 
compared to the total PYTHIA MC distributions as well as the contributions from 
the beauty, charm and light-flavour MC subsamples. All samples were normalised 
according to the scaling factors obtained from the fit (see Section^). 
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Figure 2: Distributions of (a) rf et and (b) p^ 1 of the jets associated with a sec- 
ondary vertex, (c) m vtx and (d) n tr k of the selected secondary vertices, (e) y 2 /ndf 
of the secondary vertices before the cut shown in the figure had been applied, (f) 
shows weighted by the number of jets with associated secondary vertices in the 
event. The data are compared to the total MC distributions as well as the contri- 
butions from the beauty and charm MC subsamples. All samples were normalised 
according to the scaling factors obtained from the fit (see Section^. 
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Figure 3: Distribution of the subtracted decay-length significance in three mass 
bins. The data are compared to the total PYTHIA MC distribution as well as the 
contributions from the beauty, charm and light-flavour MC subsamples. All samples 
were normalised according to the scaling factors obtained from the fit. 
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Figure 4: Distributions of (a) rf et , (b) pj, , (c) m vtx and (d) n tr k of the selected sec- 
ondary vertices and (e) subtracted decay-length significance, for a beauty- enriched 
subsample with 2 < m vtx < 7.5 GeV and \S\ > 8. The data are compared to the 
total MC distributions as well as the contributions from the beauty and charm MC 
subsamples. The light-flavour contribution is not shown separately as it is negligible 
on the scales shown. All samples were normalised according to the scaling factors 
obtained from the fit. 
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Figure 5: Distributions of (a) rf et , (b) , (c) m vtx and (d) n tr k of the selected sec- 
ondary vertices and (e) subtracted decay-length significance, for a charm- enriched 
subsample with 0.8 < m vtx < 2 GeV . No additional significance cut was applied 
here. The data are compared to the total MC distributions as well as the contri- 
butions from the beauty and charm MC subsamples. The light-flavour contribution 
is not shown separately as it is negligible on the scales shown. All samples were 
normalised according to the scaling factors obtained from the fit. 
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Figure 6: Differential beauty-jet and charm-jet photoproduction cross sections 
as defined in Section^ as a function of (a)-(b) p 3 ^ and (c)-(d) rf et . The data are 
shown as points. The inner error bars are the statistical uncertainties, while the 
outer error bars show the statistical and systematic uncertainties added in quadra- 
ture. The band represents the NLO QCD prediction, corrected for hadronisation 
effects, using CTEQ6.6 as proton PDF; the shaded band shows the estimated un- 
certainty. The NLO QCD prediction using CTEQ5M as proton PDF is depicted 
separately (dotted- dashed line). The scaled PYTHIA MC prediction (dashed line) 
is also shown. 
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Figure 7: Differential beauty-jet photoproduction cross sections as a function of 
rf et . The filled circles show the results from this analysis (the same data as shown 
in Fig. W(c)); the open circles show the results from a previously published mea- 
surement JE/. The inner error bars are the statistical uncertainties, while the outer 
error bars show the statistical and systematic uncertainties added in quadrature. 
The scaled PYTHIA MC prediction is also shown (dashed line). 
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Figure 8: (a) Summary of differential cross sections for b-quark production as 
a function of p\ as measured by the ZEUS collaboration. The measurements are 
shown as points, with the results of this analysis shown as inverted triangles. The 
inner error bars are the statistical uncertainties, while the outer error bars show the 
statistical and systematic uncertainties added in quadrature. The band represents 
the NLO QCD prediction and its theoretical uncertainty. The solid line shows 
the prediction for fi 2 = (m 2 + p^)/A, while the dashed line shows the prediction for 
fi 2 = ml+p 2 ^. (b) The ratio of the measured cross sections, a meas , to the theoretical 
prediction, o , for fi 2 = (m 2 +p^)/4. 
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Figure 9: (a) Summary of differential cross sections for c- quark jet production as a 
function ofp^ 3et (a) as measured by the ZEUS collaboration. The measurements are 
shown as points, with the results of this analysis shown as inverted triangles. The 
inner error bars are the statistical uncertainties, while the outer error bars show the 
statistical and systematic uncertainties added in quadrature. The band represents 
the NLO QCD prediction and its theoretical uncertainty. The solid line shows 
the prediction for fi 2 = (m 2 c + p|,)/4, while the dashed line shows the prediction for 
/i 2 = ml+p 2 ^. (b) The ratio of the measured cross sections, a meas , to the theoretical 
prediction, a , for ji 2 = (m 2 + p 2 -.)/4. 
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